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^ . ABSTRACT 

^ ' We investigate the photometric redshift accuracy achievable with the ^JT^i?/ infrared 

data in deep mufti-band surveys, such as in the North Echptic Pole field. We demon- 
strate that the passage of redshifted policyclic aromatic hydrocarbons and silicate 
^ ' features into the mid-infrared wavelength window covered by AKARI is a valuable 

means to recover the redshifts of star burst galaxies. To this end we have collected a 
• sample of '^60 galaxies drawn from the GOODS-North Field with spectroscopic red- 

shift 0.5 ^ Zspcc ^1 1-5 and photometry from 3.6 to 24 /im, provided by the Spitzer, 

ISO and AKARI satellites. The infrared spectra are fitted using synthetic galaxy 
Spectral Energy Distributions which account for starburst and active nuclei emission. 
For ~ 90% of the sources in our sample the redshift is recovered with an accuracy 
' l-Zphot — ^specl/ (1 + Zspec) ~ 10%. A similar analysis performed on different sets of sim- 

ulated spectra shows that the AKARI infrared data alone can provide photometric 
redshifts accurate to Izphot — 2spcc|/(l + Zspoc) ^ 10% (1 cr) at 2 ;$ 2. At higher redshifts 

the PAH features are shifted outside the wavelength range covered by AKARI and the 
, photo-z estimates rely on the less prominent 1.6 /Ltm stellar bump; the accuracy achiev- 

able in this case on (1 + z) is ^ 10 — 15%, provided that the AGN contribution to the 
infrared emission is subdominant. Our technique is no more prone to redshift aliasing 
than optical-uv photo-z, and it may be possible to reduce this aliasing further with 
the addition of submillimetre and/or radio data. 
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1 INTRODUCTION 

The discovery that infrared luminous starbursting galaxies 
are significant and possibly dominant contributors to the 
cosmic starformation history of the Universe has had an 
enormous impact on the understanding of galaxy evolution 
(Hughes et al. 1998; Lagache, Puget & Dole 2005). However, 
these galaxies are often too faint in the optical for large op- 
tical redshift compaigns, or have ambiguous optical identi- 
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fications (Chapman et al. 2004, 2005). Attention therefore 
have been focusing in the last years on infrared photometric 
redshift estimators. Sawicki (2002) showed that the 1.6 ^m 
spectral feature arising from the photospheric emission from 
evolved stars could be used to obtain crude photometric 
redshift constraints with the 3.6-8.0 fim photometry from 
the IRAC instrument on Spitzer. Unfortunately, the 8-24 /im 
wavelength gap in Spitzer's wide-field survey capability pre- 
vented the use of longer wavelength rest-frame features as 
redshift indicators, such as the much more prominent Poli- 
cyclic Aromatic Hydrocarbons (PAH) and ~ 10 /^m silicate 
absorption features. 
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AKARI is a Japanese-led infrared space telescope, 
which was launched successfully in Feb 2006. It has un- 
dertaken deep surveys near the Ecliptic Poles, far deeper 
than its all-sky survey. The North Echptic Pole (NEP) field 
is the largest deep-field legacy survey of AKARI, covering 
0.38 deg'^, and is its major deep-field legacy (Wada et al. 
2008). A wider and shallower survey at the NEP has also 
been performed with AKARI, over a 5.8 deg^ area (Lee et 
al. 2008), and further deep photometric surveys have been 
made of well-studied Spitzer fields (Pearson et al. 2008) . 

AKARI observed in 9 bands from ~2 up to ~24 fim. 
The mid-infrared imaging spans the 8 — 24 fim gap between 
Spitzer's IRAC and MIPS instruments. The unique diagnos- 
tic power of AKARFs Spectral Energy Distribution (SED) 
analysis has been confirmed by Takagi et al. (2007), who 
very successfully traced several starburst PAH features and 
identified the AGN dust torus excess over the starburst emis- 
sion, obtaining at the same time a good match between the 
photometric redshift derived from HyperZ (Bolzonella et al. 
2000) and the one obtained by fitting the infrared data with 
reliable starburst models. 

In the present work we investigate the photometric red- 
shift accuracy achievable with AKARI in deep multi-band 
surveys from infrared data alone. Infrared spectra are fit- 
ted using reference SED templates derived from the star- 
burst model of Takagi, Arimoto & Hanami (2003, hereafter 
TAII03; see also Takagi, Hanami & Arimoto 2004, hereafter 
THA04), and already exploited in Takagi et al. (2007). How- 
ever in order to deal with a possible excess over starburst 
emission due to an AGN activity we have included a set 
of reference AGN spectra derived from the model of Ef- 
stathiou & Rowan- Robinson (1995, hereafter ER95). Our 
photometric redshift code combines the starburst and AGN 
components to provide the best fit to observed or simulated 
spectra. We fit the infrared SEDs of a sample of galaxies 
with spectroscopic redshift drawn from the GOODS-North 
Field where deep AKARI observation have been performed 
at 11 and 18 fim, and of a set of simulated spectra in the 
NEP Deep Field with AKARI. We aim to demonstrate in 
this way the power of the PAH and silicate absorption feau- 
res to obtain reliable redshift estimates based on infrared 
photometry alone. 

We start in Section[2]describing the reference SED tem- 
plates. Section [3] presents the sample of infrared sources 
drawn from the GOODS-N field and Section |4] provides our 
results on the comparison between the photometric and the 
spectroscopic redshifts. Section [S] describes the simulations 
and the results on the photometric redshift accuracy for the 
AKARI 'NEP Deep field. Conclusions are summarized in Sec- 
tion [B] 



2 INFRARED SED TEMPLATES 
2.1 Starburst component 

We adopt the StarBUrst with Radiative Transfer (SBURT) 
model of THA03 and TAH04 to deal with the infrared emis- 
sion due to starburst activity. 

The model deals with the star formation and chemi- 
cal evolution in the starburst regions using the infall model 
of chemical evolution of Arimoto et al. (1992). Under the 



assumption that the amount of gas in the outflow from the 
starburst region has a negligible impact on the chemical evo- 
lution as a whole, the starburst is characterized by the rate 
of gas infall and that of the star-formation. However, for 
simplicity, THA03 set the infall time-scale equal to star- 
formation time-scale, to. The equations describing the time 
variation of gas mass, total stellar mass and gas metallicity 
are solved numerically using in input the unobscured stel- 
lar SEDs derived from the evolutionary population synthesis 
code of Kodama & Arimoto (1997). A top-heavy IMF is as- 
sumed with a slope x = 1.10, flatter than the Salpeter IMF 
{x = 1.35). The initial metallicity of the gas cloud, Zi, is 
set to zero. According to TAH03, the chemical evolution as 
a function of t/to as well as the properties of the SED are 
almost independent of the practical choice of to, speciflcally 
when to 50 Myr. Therefore, all the reference SED tem- 
plates used here for the starburst are calculated adopting 
to = 100 Myr. 

THA03 assumes a starburst region in which stars and 
dust are distributed within a radius -Rt, and introduce the 
following mass-radius relation 

^^ef^I^Y (1) 

Ikpc l^lOSMey ' ^' 

where O is a compactness factor that expresss the mat- 
ter concentration, with the mean density being higher for 
smaller values of O, while M* is the (time-dependent) to- 
tal stellar mass in the starburst region. The exponent 7 is 
set to ^ which results in a constant surface brightness of 
the straburst galaxy for constant 0. Dust is assumed to be 
homogeneously distributed within Rt while a King profile 
is adopted for the stellar density distribution. The optical 
depth of the starburst regions is a function of the starburst 
age and of the compactness factor, and increases for decreas- 
ing values of O. 

The amount of dust is given by the chemical evolution 
assuming a constant dust-to-metal ratio. So , with three dif- 
ferent dust models, i.e. the model for dust in the Milky 
Way (MW), the Large MageUanic Cloud (LMC) and the 
SmaU Magellanic Cloud (SMC). The values of So derived 
from the extinction curve and the spectra of cirrus emission 
are 0.40, 0.55 and 0.75 for the MW, LMC and SMC, respec- 
tively. The difference among the three extinction curves is 
attributed to the variation of the ratio of carbonaceous dust 
(graphite and PAHs) to silicate grain with the fraction of 
silicate dust grains increasing from MW to SMC model, i.e. 
as a function of metallicity; therefore the silicate absorption 
features are most prominent in the SEDs described by the 
SMC type dust model. The SED from ultraviolet to submil- 
limetre wavelengths of starburst region is calculated for each 
starburst age, compactness parameter and extinction curve 
using a radiative transfer code that assumes isotropic multi- 
ple scattering and accounts for self-absorption of re-emitted 
energy from dust. 

The absorption/emission properties of the dust are re- 
sponsible for the infrared SED features of starforming galax- 
ies: the PAHs emissions at 3.3, 6.2, 7.7, 8.6 and 11.3 /xm 
and the silicate resonances at 9.7 and 18 ^m. An example 
of starburst SED template drawn from the Takagi 
et al. model, shoviring the main infrared features, is 
presented in Fig. [TJ The efficiency of these features as 
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Figure 1. Example of a starburst SED template (rest-frame) 
drawn from the Takagi et al. model, showing the main infrared 
features exploited as redshift indicators: PAH emissions at 3.3, 
6.2, and 7.7 ^m, and silicate absorption at 9.7/im. The 1.6 ^m 
stellar bump is also visible. In the central and in the bottom 
panels are shown the transmission curves of the infrared filters 
exploited in the GOODS-N field and in the AKARI NEP deep 
field, respectively, normalized to a majdmum value of one: Spitzer 
3.6, 4.5, 5.8, 8.0, 16 and 24 /im filters ( green), ISO 6.5 and 15 ^tm 
filters (blue), AKARI Slters (the 11 and IS/im filters in yellow, 
the others in light blue). 



redshift indicators is what we are going to put to the test in 
the present work. 



2.2 AGN component 

The infrared spectrum of the AGN dusty torus is modelled 
as in ER95. They combine an accurate solution for the ax- 
ially symmetric radiative transfer problem in dust clouds 
with the multigrain dust model of Rowan- Robinson (1992) 
and three different geometries for the dusty torus that sur- 
round the central supermassive black hole. Here we assume 
their thick tapered disc model, which is a disc-like torus 
whose height increases with the distance from the central 
source but tapers to a constant height in the outer parts. 
This choice has been found to provide the best agreement 
with the observational constraints on AGN. The tapered 
disc is assumed to have a density distribution. The ra- 
tio between the outer and the inner radius of the torus is set 
to 20, in the middle of the values considered by ER95, while 
we assume a value of 45° for the half-opening angle of the 
torus (the values assumed by ER95 for the nucleus of NGC 
1068). The equatorial optical depth of the torus is fixed to 
1000. Therefore the SED of the AGN torus depends only on 
the viewing angle of the torus, dview ■ 

With the above choice of parameters the SED shows deep 
absorption features at ~ 10 fim due to silicate dust for edge- 



on views of the torus {O^iew = 0°), but it is rather featureless 
when the torus is seen face-on {9 view ~ 90°). 



2.3 Reference SED templates and fitting 
procedure 

Using the models described above we have created a set 
of reference SED templates meant for fitting any set of in- 
frared photometric data (both real and simulated). These 
templates have been derived as follows (see TAH04). 

• STARBURST component: for each type of extinction 
curve, we adopt 10 different starburst ages in the interval 
t/to = 0.1 - 6.0 (or equivalently t = 10 - 600 Myr, being 
to ~ 100 Myr), and 17 different compactness factors in the 
range Q = 0.3 — 5. 

• AGN COMPONENT: we vary the viewing angle between 
0° and 90°, with 11 discrete values. 

The resulting SED templates are convolved with the re- 
sponse functions of the instruments for different values of 
the redshift in the range z = — 7 and in steps of 0.02. The 
fluxes derived for the two SED components (i.e. starburst 
and AGN torus) are then linearly combined to give the pre- 
dicted flux at the passbands of the instruments. The best flt 
parameters of the SED models, including the redshift, are 
then obtained by minimizing the 



2 

X = 



E 



(2) 



where Fahs{i) is the measured flux in the passband i and a[i) 
is its corresponding uncertainty. Fk{i) represents the pre- 
dicted flux in the passband i contributed by the SED com- 
ponent k, with fc = 1, 2 denoting "Starburst" and "AGN", 
respectively, fk is therefore the relative contribution of the k 
component to the total bolometric luminosity of the galaxy. 
Note that for each set of values of the SED model param- 
eters, the values of /i and /2 are fixed by the condition of 
minimization with the restriction that they must be all non- 
negative. When a negative value is obtained for either /i or 
/2 , then that parameter is set to zero and the process of min- 
imization is performed on the other SED component alone. 
We descard those sets of values of the SED model parame- 
ters for which the bolometric luminositjQ resulting from the 
minimization over the parameters /i and /2 hes above 10^* 
L0 or below 10* L©. Since we are dealing with photometric 
data from different telescopes we have set a minimum error 
of 10% on the measured fluxes before performing the SED 
fltting. 

The goodness of the fit is described by the probability for 
a x^-distribution, with the number of degree of freedom set 
by the problem, to produce a value of higher than the 
one obtained by the best SED-fit. The number of degree 
of freedom is given by = Ji-data — JT-par where Tidata is the 
number of fitted data (ranging from 8 to 10 depending on 
the availability of ISO and/or IRS photometry) while ripar is 
the number of parameters of the model, i.e. the redshift, the 
normalization of the SED components and the SED model 



^ The bolometric luminosity is obtained by integrating the SED 
in the rest-frame wavelength interval 0.1-1000 /.tm. 
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Table 1. Some information about the Spitzer and AKARI in- 
frared catalogues used in the present work: observing wavelength, 
Full With at Half Maximum (FWHM) of the instrument, number 
of surces detected above 3cr, corresponding minimum source flux, 
and surveyed area. 



instrument waveband FWHM tl sources S^nin Area 

{^m) (aresee) {^Jy) {arcmin^) 



Spftzer/IRAC 3.6 1.6 5792 0.52 230 

Spitzer/mAC 4.5 1.7 5576 0.53 230 

Spitzer/mAC 5.8 1.9 2328 2.74 226 

Spitzer/mAC 8.0 2.0 2186 1.80 249 

AKARI/mC 11.0 4.8 242 44.04 101 

AKARI/mC 18.0 5.7 233 96.29 115 

Spitzer/MIPS 24.0 6.0 1199 80.00 254 



the 6.5fim band. It includes 49 objects, 42 detected only at 
15 /xm, 3 at only 6.5 and 4 at both wavelengths. Aussel et 
al. also provide an additional, less secure, list of 51 sources 
of which 47 are detected at 15 fim only, 4 at 6.5 fj,m only, but 
none in both filters. All together the two catalogues include 
a total of 100 objects. 



3.2 Spitzer 

The Spitzer data for the HDF-N are part of the GOODS 
Spitzer Legacy Data Product^ and are in the public do- 
main. The Spitzer data sets include the images of both 
GOODS fields at 3.6, 4.5, 5.8, 8.0 /im from the Infra-Red Ar- 
ray Camera (IRAC) and at 24 /im from the Multiband Imag- 
ing Photometer for Spitzer (MIPS), plus a list of sources for 
the MIPS 24 ^m imaging. The MIPS catalogue consists of 
1199 sources detected at 3a, with fiux densities greater than 
80/iJy (see Table [TJ, a limit where the source extraction is 
stated to be highly complete and reliable. Source lists for the 
IRAC imaging have not yet been released by the GOODS 
consortium. Therefore we derived our own IRAC source cat- 
alogues in GOODS-N using SExtractor (Berlin & Arnouts 
1996) with the default set of values for the configuration 
parameters. The source fiuxes and the corresponding errors 
have been derived from the automatic aperture magnitudes. 
Sources were extracted down to a minimum signal-to-noise 
ratio of 3. The resulting catalogues consist of 5792, 5576, 
2328, 2186 sources at 3.6, 4.5, 5.8 and 8.0 fim, respectively 
(see Table g. 

Imaging at 16 fim of selected areas within the GOODS-N 
field have been obtained as a result of a pilot study with 
the Spitzer IRS (Teplitz et al. 2005). The majority of the 
area (30 of 35 arcmin^) reaches a 3o" sensitivity of ~75/iJy. 
Teplitz et al. (2005) detected 149 objects at 16 fj,m, with 
flux ranging from 21 jj,Jy to 1.24 mjy, and with photometry 
in good agreement with the 15 /im ISO survey of the same 
area. 



3.3 AKARI 

The GOODS-N region has recently been imaged at 11 ^m 
and 18 /im with the AKARI Infra-Red Camera (IRC) within 
the FUHYU program. FUYHU is an AKARI mission pro- 
gram to follow-up well-studied Spitzer fields (Pearson et 
al. 2008) including the GOODS-N field. This survey has 
mapped the Lockman-Hole and ELAIS Nl fields, other well- 
studied fields with sufficiently high ecliptic latitudes and 
consequently high AKARI visibility. The source extraction 
and fiux calibration are described in Pearson et al. (2008), 
and are based partly on the source extraction methodology 
developed for sub-millimeter surveys (Serjeant et al. 2003). 
The samples comprises 242 detections at ll^m and 233 at 
18 /xm respectively, with a signal-to- noise ratio higher than 
3 (see Table [TJ. 



^ see http:/ /data. spitzer. caltech.edu/popular/goods/Documents/ 
goodsdataproducts.html 



parameters. As a resul10 ripar = 6 if both /i and /2 are pos- 
itive, otherwise ripar = 4 (/2 = 0) or Tipar = 3 (/i = 0). A 
99% confidence interval on the estimated photometric red- 
shift is derived by the Ax^ method (Avni 1976), being the 
number of "interesting parameters" equal to 1, i.e. Zphot- 
In this case the value of Ax^ defining the confidence inter- 
val at 99% is 6.63. A best SED-fit is considered "good" if 
P^2 > 1%. 



3 DATA SETS 

In order to show the reliability of the redshift estimates 
based on the PAH and silicate infrared features, we have 
assembled a sample of galaxies with flux measurements in 
the wavelength range 3.6-24 fim and with measured spec- 
troscopic redshifts. The sample has been selected within the 
northern fleld of the Great Observatories Origins Deep Sur- 
vey (GOODS, Dickinson et al. 2001) because of its richness 
in multiwavelength photometric data and spectroscopic red- 
shifts. The GOODS observations covered two flelds on the 
sky: a northern target area (GOODS-N) coincident with 
(but significantly larger than) the Hubble-Deep Field North 
(HDF-N, Williams et al. 1996), and a similarly sized south- 
ern field (GOODS-S) coincident with the Chandra Deep 
Field South (Giacconi et al. 2000). The GOODS-N field 
has been imaged at infrared wavelengths by different in- 
struments: at 6.75 //m and 15 /xm by ISOCAM on board of 
the Infrared Space Observatory (ISO, Kessler et al. 1996), 
at 3.5, 4.5, 5.8, 8 fj, and 24 fira by Spitzer, and, more re- 
cently, at 11 and 18 ^im by AKARI (Pearson et al. 2008). 
Some subregions of the GOODS-N field (including part of 
the HDF-N) were also imaged at 16 ^m by the Spiteer Infra- 
Red Spectrograph (IRS) blue peak-up filter (Teplitz et al. 
2005). It is these sets of infrared data we have exploited here 
and we provide brief descriptions below. 

3.1 ISO 

We use the catalogue of ISO sources in the HDF-N produced 
by Aussel et al. (1999). The source catalogue is claimed to be 
95% complete at 200 /iJy in the 15 fim band and at 65 /xJy in 

2 The type of the extinction curve (i.e. MW, LMC, SMC) is not 
considered as a SED parameter and therefore it has not been 
included in the calculation of 
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3.4 Spectroscopic redshifts 

The region around tiie HDF-N liad been tlie subject of in- 
tensive spectroscopic campaigns in tlie 90s' by a variety 
of groups using the Low Resolution Imaging Spectrograph 
(LRIS; Oke et al. 1995) on the Keck 10 meter telescopes. A 
compilation of the LRIS-spectra for 671 sources is presented 
in Cohen et al. (2000). Subsequently, two parallel spectro- 
scopic projects were carried out in the GOODS area using 
the Deep Imaging Muhi-Object Spectrograph (DEIMOS) 
on the Keck II telescope. The Team Keck Treasury Sur- 
vey (TKTS) of the GOODS-N field (Wirth et al. 2004) 
has focused on an Rab < 24.4 magnitude selected sam- 
ple while the survey by Cowie et al. (2004) was embedded 
within observations targeted at high-redshift galaxies and X- 
ray and radio selected galaxies. More recently, Steidel and 
collaborators have conducted a spectroscopic survey in the 
GOODS-N with the blue arm of LRIS (LRIS-B) targeting 
several hundreds of star-forming galaxies and AGNs at red- 
shifts 1.4 ^ z ^ 3.0. The source candidates for spectroscopic 
follow-up were pre-selected using different criteria based on 
their optical colours (Steidel et al. 2003, 2004; Adelberger et 
al. 2004). The resulting spectroscopic catalogue, presented 
by Reddy et al. (2006) includes 342 objects and provide for 
each of them both optical photometry and infrared photom- 
etry at 3.6, 4.5, 5.8, 8.0 and 24 /^m from Spitzer. We use here 
the compilatiorlfl of redshifts assembled by the Keck Team 
combining all the spectroscopic surveys within the GOODS- 
N field produced up to the 2004, and we add to it the 
spectroscopic sample of Reddy et al. (2006). Note that the 
TKTS is based primarly on observations with the DEIMOS 
instrument on Keck II, whose spectral range only allows 
the detection of emission and absorption features from ob- 
jects at 2 ^ 1.2. Conversely, the selection criteria used in 
the spectroscopic survey of Reddy et al. is better at identi- 
fying galaxies at z 1.4. We found indeed that only 59 out 
of the 342 spectroscopic sources listed in the Reddy et al. 
catalogue are already included in the TKTS compendium 
redshift catalogue. 



bands over the whole field and in HK' over a smaller re- 
gion covering the Chandra Deep Field South, down to 5a 
(AB magnitudes) limits of 27.1, 26.9, 26.8, 26.6, 25.6, 25.4 
and 22.1, respectively. The images and the corresponding 
catalogues are available on the World Wide Wellf|. 

Here, optical data are used just for comparison with the 
infrared imaging, but they are not exploited in the SED fit. 



3.5 Optical data 

Optical imaging of the GOODS-N field was obtained with 
the Avanced Camera for Survey (ACS) on board the Hubble 
Space Telescope (HST) at B, V, i', z' bands. The ACS images 
and the source catalogues extracted by the GOODS Team 
are of public domairlf|. 

Optical/near-infrared ground-based images are also 
available for the GOODS-N field (Capak et al. 2004). An 
intensive multi-color imaging survey of ~ 0.2 deg'^ centered 
on the HDF-N have been carried out using different instru- 
ments: the Kitt Peak National Observatory (KPNO) 4 meter 
telescope with the MOSAIC prime focus camera, the Subaru 
8.2 meter telescope with the Suprime-Cam instrument, and 
the QUIRC camera on the Hawaii 2.2 meter telescope. The 
surveyed area is referred to as the Hawaii HDF-N (Capak 
et al. 2004). Data were collected in U, B, V, R, I and z' 



* http:/ /www2. keck. hawaii.edu/tksurvey/ 
^ http:/ /archive. stsci.edu/prepds/goods/ 



http:/ /www. ifa.hawaii.edu/ capak/hdf/index.html 



Table 2. List of infrared sources with spectroscopic redshift drawn from the GOODS-N field, by cross-matching the source lists of Aussel et al. (1999) and of Teplitz et al. (2005) with 
the Spitzer and AKARI catalogues in the same field. The source positions reported here are those of the spectroscopic counterpart. 
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Table 2. Continued 
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3.6 Cross-matching 

We have cross-matched the spectroscopic catalogues with 
each of the infrared catalogues mentioned above by using 
for the matching radius a value equal to twice the Gaussian 
rms width apwHM (=FWHM/2V2 In 2) of the instrument 
beam, where ufwhm is 6.4" for ISOCAM at 15/im, 0.68", 
0.72", 0.81", 0.85" for IRAC at 3.6, 4.5, 5.8, S.O^im respec- 
tively, 2.55" for MIPS at 24^im, 1.53" for IRS at 16 ^im, 
2.04" and 2.42" for the AKARI at 11 and 18^im, respec- 
tively. In order to ensure full spectral coverage from mid- to 
far-IR wavelengths we kept only spectroscopic sources with 
a counterpart in all of the Spitzer and AKARI mfTaxed cat- 
alogues considered here plus a counterpart in either the ISO 
15 /Km or the IRS 16 fim catalogues. The advantage of cov- 
ering the whole mid- to far-IR spectral range is made clear 
by looking at the middle panel of Fig. [T] where the Spitzer, 
AKARI and ISO filters are shown: all the main PAH and sil- 
icate features are sampled up to 2: ~ 2 with the photometric 
data exploited here. At z ^ 2 the PAH features are shifted 
outside the wavelength range covered by Spitzer and AKARI 
so that they can no longer be exploited as redshift indica- 
tors. In this case, the rest-frame 1.6 ^m bump can be used 
instead for photometric redshift estimates (Sawicki 2002). 

We ended up with a sample of 59 objects; 22 of them 
have flux measurements at both 15 /im and 16 fim. ISO 
6.5 fim flux measurements (or upper limits) have been in- 
cluded, when available, in the SED fitting process. The red- 
shifts of the sources and their infrared fluxes are listed in 
Table [2I We also provide in Fig. [2] postage stamp images for 
all the sources in our sample at R, HK' , Spitzer and AKARI 
wavebands. 

For almost all of the sources in the sample we found a 
single counterpart at each infrared waveband. However few 
of them appear to lie in crowded regions, possibly affect- 
ing the infrared flux at the longest wavebands, see e.g. ID3, 
ID20, ID25, ID31, ID35, ID37, ID38 and ID44. There are in 
particular pairs of objects "sharing" the same AKARI (and 
ISO) fluxes, i.e. ID31-ID35 and ID37-ID38. In this cases the 
AKARI (and ISO) photometry should be taken as an upper 
limit to the intrinsic flux of the source at those wavelengths. 

Most of the sources in the sample lie at z ^ 1.5. Only 

3 out of 59 have ^^poc > 1.5: ID5 (z^pcc = 2.5920), ID41 
(zspec = 2.0050) and ID59 (z^pec = 2.2110). We found that 
just five objects in the Reddy et al. catalogue have a coun- 
terpart at both 11 and 18 /xm and either at 15 or 16 fim. 
They are (according to the names they have in the Reddy 
et al. source list): BX1321 (at z=0.139), BM1156 (at 2.211), 
BM1299 (at z=1.595), BM1326 (at 1.268), and MD39 (at 
2.583). Four of them are in common with the TKTS spectro- 
scopic catalogue and fall into the final sample of 59 infrared 
sources presented here, corresponding to IDl, ID5, ID52 and 
ID59. MD39 is instead missing from our final sample because 
it lacks a counterpart in our IRAC catalogues at both 3.6 
and 8.0 fim. MD39 appears to be the source very close to 
ID24 in the corresponding postage stamp image of Fig. (2] at 
wavelength A ^ 5.8 fim. The pair of objects is not resolved 
at longer wavelengths and indeed Reddy et al. do not pro- 
vide a 24 fim flux estimate for it, probably because it cannot 
easily be deblended from the nearby ID24 source (although 
they provide for it flux estimates at all the four IRAC wave- 



bands). For this reason we have decided not to include MD39 
in our final catalogue and to add instead ID24 to the list of 
our infrared sources with uncertain far-IR photometry. 

Fig. |3] shows the observed flux ratios as a function of 
the spectroscopic redshifts (black dots), and compare them 
to the predictions from the starburst reference templates 
for a representative set of SED model parameters (blue 
curves), and from the AGN reference templates (megenta 
curves). The effects of the PAH 6.2 nm emission feature and 
of the 9.7 /im silicate absorption are clearly manifest in the 
data. The deep "well" observed in the Fs.g/Fg.o, Fg.o/Fn, 
F11/F15 and F11/F16 diagrams at 2; ~ 0.2, z ~ 0.6, z ~ 1.2 
and z; ~ 1.4 is due to the passage of the PAH 6.2 fim feature 
through the 8.0, 11, 15 and 16 nm wavebands, respectively. 
The 9.7 fim silicate absorption manifests itself as a bump in 
the 8.0/11 fim flux ratios around z ~ where the feature 
falls into the AKARI 11 /im passband. As the same feature 
enters the 15 and 16 /^m wavebands, which occurs at z ~ 0.5, 
it produces a significant corresponding peak in the F11/F15 
and F11/F16 diagrams. A hint of a bump due to the silicate 
absorption is seen also in the 15/18 /im and 16/18 fim fiux 
ratios at 2 ~ 0.7, although it is not particularly prominent 
due to the wider wavelength coverage of the AKARI 18 fim 
filter compared to the 11 fim band. Finally, by entering the 
MIPS 24 fj,m band, the silicate absorption induces another 
significant bump in the 18/24 fim flux ratio at z ~ 1.5. Note 
that such a bump can be used for efficiently selecting ultra- 
luminous infrared galaxies in the redshift range 1 ^ z ^ 2 

with AKARI (see Takagi & Pearson 2005). The ~ 10 ^im sil- 
icate absorption is also responsible for the behaviour in the 
predicted flux ratios of the AGN templates when the torus 
is seen edge-on. For a face-on AGN the corresponding flux 
ratios are almost independent of redshift, implying that for 
a power-law infrared spectrum the recovery of the redshift 
from mid-/far-IR photometry alone is extremely challeng- 
ing, if not impossible. 

On average the range of flux ratios spanned by the data 
is accounted for by the model with the sole exception of the 
15/16 fim flux ratios. Indeed we observe for few objects a sig- 
nificant steep increase of the flux from 15 to 16 /im (ID28, 
ID39, ID42, which lie well below the theoretical expecta- 
tions) or, conversely, a notable decrease of the flux when 
moving from 15 to 16 ^m (ID44). 



4 RESULTS OF SED FITTING 

The results of the SED fitting are listed in Table [3] A com- 
parison between the photometric data and the best-fit SED 
model is shown in Fig. 3] The derived photometric redshifts 
versus the spectroscopic redshifts are presented in Fig. [5] 
where sources with P^2 > 1% are indicated by filled cir- 
cles (42 in total, i.e. ~ 70% of the whole sample), and 
those with P^2 < 1% are identified by open circles. In the 
same figure the solid line marks the ideal case in which 
2phot = Zspcc, while the dot-dashed lines delimit the region 
where l^phot ~ 2apcc|/(l -I- ^spoc) < 10%. Error bars corre- 
sponding to a 99% confidence limit have been drawn, for 
reason of clarity, only for objects with P^2 > 1% which lie 
within the 10% confidence region. 

For almost all of the sources with Zspcc ^ 1.5 and 

P^2 > 1% the redshift is recovered with an accuracy 
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Figure 2. 20" X 20" postage stamp images centered on the infrared sources listed in Table [2] From left to right: R and HK' optical 
images, IRAC images at 3.6, 4.5, 5.8 8.0 /xm, /lifylij/ images at 11 and 18 /^m, and MIPS image at 24 /im. Red dots in the Spitzer a.Tid 
yli<'j4ij/ images represent 3(7 detections at the corresponding wavebands while the circle is the "error" -circle of 2o'pwHM-radius used for 
cross-matching the spectroscopic catalogue with the infrared catalogues, i.e. 2o'pwHM = 4.08", 4.84" and 5.10" at 11, 18 and 24 ^m, 
respectively. A circle of 1.5" in the R band image indicates the spectroscopic counterpart. 
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Figure 2. Continued. 
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Figure 3. Observed flux ratios as a function of spectroscopic redshift (black dots) for the sample of infrared sources extracted from the 
GOODS-N field. Blue lines show the predicted colors from the starburst models of Takagi et al. for a subsample of representative values 
of the sed model parameters (see text). The predicted colors from the AGN templates exploited in the SED fitting are shown in magenta. 



I^phot — 2spcc|/(l + ^spcc) ^ 10%. The only exceptions are 
ID22 (at Zspoc = 0.5560), ID42 (at z^pec = 0.5480) and ID44 
(at Zspec = 1.4865). 

For ID22, we would expect the PAH 6.2 and 7.7 ^im 
features to enter the 11 /.im AKARI filter at Zspoc ~ 0.5 and 
therefore to produce a bump in the measured flux at that 
waveband. This is for example what we observe in the spec- 
trum of ID7 and ID 15 whose redshifts are very close to that 
of ID22: a bump at ll^m followed by a plateau at longer 
wavelengths as a result of the convolution of the instruments 
filters whith the redshifted silicate absorption features at 9.7 
and 18 /im and with the PAH 11.3 /im mission feature. On 
the contrary, the SED of ID22 increases monotonically in 
the interval 8-18 /xm. Interestingly, we found that at 2 ~ 0.5 
the highest contribution to the value of the in Eq. [2] 
(after minimization over the other SED model parameters) 
comes from the AKARI 11 /.im waveband. By removing the 
data point at 11 iim from the measured spectrum and per- 
forming the fit on the other photometric data alone we get 
Zphot = 0.56, in excellent agreement with the spectroscopic 
value. We also tested that none of the flux measurements 
at wavelengths > 11 ^ra have such an effect on the photo-z 
estimate. Therefore the failure in the recovery of the redshift 
for ID22 seems to be determined exclusively by the AKARI 
photometry at 11 /im. This result points to the conclusion 
that the source has an intrinsically low PAH 7.7 /im emission 
that our reference SED templates are not able to account for 
that. 

ID42 is almost at the same redshift of ID22 and its spec- 
trum has a peculiar power-law shape above 8 /im. The ISO 
data suggest an emission feature at 6.5 ^m that, if due to 
the restframe 3.3 /im PAH feature, would suggest a redshift 
2: ~ 1 at variance with the spectroscopic value. However the 
measurement is affected by large uncertainties. By fitting the 
measured SED with the redshift fixed to its spectroscopic 
value we get a minimum reduced of ~ 3.5. Examining 
the contribution of the individual wavebands to the value of 



the minimum we find the AKARI and the 4.5 /im Spitzer 
photometry to have the highest discrepancy between obser- 
vation and theoretical predictions. The fluxes at 4.5, 11 and 
18 appear lower than expected. These results indicate the 
limitation of our reference SED templates in accounting for 
the infrared spectrum of ID42. Therefore ID42 represents a 
challenge for the SED models exploited here. 

Failure in recovering the redshift of ID44 is quite 
certainly due to flux "contamination". In fact, this object 
lies in a crowded field (see Fig. [Sjl: there are two sources 
individually detected by Spitzer (with comparable fluxes at 
24 /xm) which are blended at AKARI (and ISO) wavebands. 
The spectrum of ID44 displays two bumps at 15 and 18 /im 
that the SED fltting interprets as redshifted PAH 7.7 and 
11.3 /.im features, respectively, thus implying Zphot < 1- 
Adopting the flux measurements at 15 and 18/.im as upper 
limits does not improve the flt because the SED would 
assume a featureless power-law shape. Therefore the only 
way to obtain a better photo-z estimate for this object is 
to use optical and radio photometry with higher spatial 
resolution. 

For sources at Zspcc ~ 2 the photometric redshift is 
found to be systematically and signiflcantly lower than the 
spectroscopic value, i.e. Zphot ^ 0.5. This is not surprising 
if examine the measured spectrum of the three objects 
(see Fig. The SEDs of ID5 and ID41 are completely 
featureless, resembling a power-law functional form, thus 
making the estimate of the redshift impossible without 
the support of photometric data at other wavelengths. 
On the other hand the SED of ID59 seems to indicate 
an absorption at ~ 15 /im (apart from that the SED is 
very close to a power- law). Our SED-fitting procedure 
interprets this feature as due to the redshifted 9.7 /xm 
silicate absorption in the spectrum of a young starburst, 
providing Zphot ~ 0.5. In fact, at the spectroscopic redshift 
of the source (i.e. Zspcc ~ 2.2), the 9.7 /^m absorption 
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Figure 4. Comparison between the best-fit SED model (black solid curve) and the measured infrared spectrum (dots with error bars) for 
the infrared sources listed in Table [2] The inferred contributions from a starburst and an AGN torus, are represented by the dashed and 
the dot-dashed curves, respectively. Also shown are the spectroscopic redshift of the sources (zepec), the derived photometric redshifts 
(■2phot)i the corresponding value of the minimum reduced ix^in/'^) associated probability (P^2). The red solid curve is the 

best-fit SED model obtained by assuming a priori that the AGN component is null. 
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Figure 4. Continued. 
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Table 3. Best SED-fit parameters obtained for the sample of infrared sources listed in Table[2] photometric redshif (2photi second column) 
with errors corresponding to a 99% confidence limit, minimum reduced ix^ ■ /^i third column), number of degrees of freedom {y, 
fourth column), probability associated to the minimum (^x^' ^^^^ column), age of the starburst (Age, sixth column), compactness 
factor (0, seventh column), extinction curve (Ext., eighth column), viewing angle of the AGN dust torus (^viewi ninth column), luminosity 
contributed by the starburst component (logLsBi tenth column), luminosity contributed by the AGN component (logZ/AGNi eleventh 
column). The best- fit luminosities have been derived integrating the rest- frame SED from 0.1 to 1000 lira. The spectroscopic redshifts of 
the sources (2:spec) are listed for comparison. The other columns show the results obtained by setting o, pviovi the AGN fraction, f2^ to 
zero and fitting the observed SED with the starburst SED templates alone. An asterisk marks the sources with uncertain counterpart at 
AKARI and ISO wavebands. 
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Figure 5. Photometric versus spectroscopic redshift for the sources listed in Table|2] Filled circles mark the cases of a "good" fit to the 
SED (i.e. P^2 > 1%) while the open circles correspond to a bad SED fit (i.e. P^2 < 1%). Sources with uncertain photometry at AKARI 
and ISO wavebands are circled in color. The dot-dashed lines delimit the region where (^phot ~ ^spec)/(l + z^spec) < 10% while the solid 
line corresponds to the case 2phot = -^spec- Errors on photometric redshifts are in 99% confidence limit and, for reason of clarity, they 
are shown only for the objects with P^2 > 1%, lying within the 10%-accuracy region. 



falls above the 24 waveband covered by MIPS. It is 
therefore plausible that the observed absorption is due 
to the transition between the rest-frame 3.3 and 6.2fim 
PAH features. In that case the source would have a much 
higher intrinsic luminosity than the value we obtained for 
-Zphot = 0.54. An independent constraint on the infrared 
luminosity of the source, e.g. from measurements at sub- 
millimeter and/or decimetric radio wavelengths, would be 
crucial for better recovering the redshift of this source. 
We plan to test the addition of extra data on the photo-z 
estimates (and on SED model parameters) in a future paper. 

We note that the estimated 99% error on the photo- 
metric redshift is particularly large for few sources lying in 



the redshift range 0.8 < «apec ^ 1, i.e. ID32, ID43, ID48, 

and ID56. This follows from the shape of the SED of these 
objects. For ID32, ID43 and ID48 the measured spectrum 
is close to a power law for A ^ 8/im, and this induces a 

degeneracy between the redshift and the normalization of 
the AGN sed component (/2) around the minimum value of 
the x^. The case of ID56 is more peculiar. At the redshift 
of the source (i.e. z^pcc ~ 1.017), the rest-frame 9.7 /^m 
silicate absorption enters the AKARI 18 ^im filters, which 
is the broadest among the AKARI mters (FWHM~10 //m); 
as a consequence the measured SED appears relatively 
flat between 15 and 24 /im. However, while for objects like 
ID21, ID30, ID36, ID45, ID46 (all at 2 = 1), a hint of the 
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Table 4. 5a flux limits for the AKARI NEP-Deep Survey used 
in the simulations (from Wada et al. 2008). 
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Figure 6. Comparison of the photometric redshift accuracy 
achieved by fitting the infrared spectrum of the sources listed 
in Table[2]with 2 SED components (i.e. starburst + AGN; upper 
panel) and with the starburst templates alone (lower panel). The 
meaning of the lines is the same as in Fig. \5\ 



silicate absorption at 18 is still visible, in the case of 
ID56 the SED is featureless above 15 fim. Consequently, the 
observed spectrum can be still reasonably well fitted by a 
mixture of starburst and AGN emission in a large redshift 
interval around z ~ 1. 

According to our SED fitting results, sources with 
P^2 > 1% and |5z|/(H-Zspec) ~ 10% manifest different levels 
of AGN emission in their SED. In same cases the SED is well 
described by a starburst template (e.g. ID17, ID35, ID38), 
while in other cases an underlying AGN activity is found to 
contribute to the measured fluxes either at all wavebands 
(e.g. ID24, ID42, ID48, ID52, ID58) or only at longest wave- 
lengths, i.e. 18 and 24 /xm (e.g. IDS, ID16). There are also 
some very peculiar/interesting cases of sources at Zspcc ~ 1 
for which the flux above A ~ 10 ^m is completely or signifi- 
cantly accounted for by the emission from an edge-on viewed 
dusty torus (e.g. ID21, ID27, ID30, ID39, ID46). 

It is worth noticing that our SED-fitting procedure is 
meant to combine the SB and AGN SED components in 
such a way to minimize the value of the in Eq. [2] and 
not to maximize the probability associated to the minimum 
value of the ■ In fact, while the value of P^2 depends on 
the number of degree of freedom of the problem (and on 
the assumption that the quantity in Eq. [2] follows exactly a 
X^-distribution), that of the minimum does not. Conse- 
quently, there may be some sources for which the best fit to 
their SED gives a smaller value for P^2 when the fit is per- 
formed by setting a priori f2 to zero than when both /i and 
/2 are taken as free parameters, even if the resulting min- 
imum x'^ is smaller in the latter case. In Table [3] we show, 
for comparison, the minimum reduced and the associated 
probability P^2 obtained when the measured SED is fitted 
with the starburst templates alone. 
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We note that there are several cases for which the value of 
P^2 is slightly lowered when using two SED components and 
in two cases it drops below the 1% limit (ID39, ID40). On the 
other hand, accounting for an AGN component improves the 
goodness of the fit for many other objects in the sample, par- 
ticularly when the measured SED is close to a power-law (see 
e.g. IDS, ID41 and ID48). However, for these extreme cases, 
a better fit doesn't imply a better redshift accuracy which 
instead remains quite poor. In Fig. [6] we compare the photo- 
metric redshift accuracy, (zphot — Zspcc)/(1 + Zspoc), achieved 
by our two-SED fitting procedure (upper panel) with that 
obtained by fitting the observed SED with the starburst 
templates alone (lower panel). In both panels filled circles 
mark the sources with P^2 > 1% while open circles indi- 
cate sources with P^2 < 1%. There are few cases in which, 
the redshift of the source is better recovered when account- 
ing for an AGN component: 1D24, 1D29, ID43, however, we 
do not see on average a significant difference in the achieved 
redshift accuracy between the two-SED fitting approach and 
that based on the starburst templates alone. 

For comparison, we also show in Fig U the result of 
the best-fit obtained when /2 — 0. There are some sources, 
mainly at redshift z < 0.5, e.g. IDl, ID2, ID4, 1D12, ID13, 
ID47 and 1D57, for which the fit to the SED is very poor 
(i.e. P^2 <^ 1%) independently on the adopted SED fitting 
approach. Failure in reproducing the infrared spectrum of 
these sources could be due to an underestimate of the fluxes 
at 11 and 18 /im. In fact, AKARI photometry were obtained 
under the assumption of point sources (see Pearson et al. 
2008 for details), while most of the objects listed above are 
significantly extended in optical images (see Fig. [2| . More- 
over, IDl and 1D47 appear to be elliptical galaxies whose 
spectrum is not expected to be well-recovered by the star- 
burst SED templates exploited here. 



5 SIMULATIONS OF THE AKARI NEP DEEP 
SURVEY 

In order to investigate the photometric redshift accuracy 
achievable with the AKARI NEP Deep Survey we have gen- 
erated three different sets of 5000 simulated spectra, includ- 
ing contributions from both starburst and AGN. 

In all the simulations, which are described in the next 
subsections, the derived starburst and AGN components are 
linearly combined together with a relative contribution to 
the bolometric luminosity randomly selected between and 
1. Redshifts are randomly assigned in the range 0-5, assum- 
ing a uniform distribution, and the spectra are then normal- 
ized to the 24 //m fluxes generated from the source count 
model of Lagache et al. (2003). A minimum 24 ^m flux of 
100 /^Jy is assumed. The spectra are then convolved with the 
nine AKARI/IRC filters, spanning the wavelength range 2- 
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Figure 7. Simulation of photometric redshift accuracy with the AKARI NEP Deep Survey. The simulated spectra have been derived 
from the same models used for constructing the reference SED templates, i.e. Takagi et al. (2003, 2004) for the starburst and Efstathiou 
& Rowan-Robinson (1995) for the AGN. The filled circles with error bars are the mean value and the Icr dispersion provided by the 
Gaussian fit to the histogram of the ^phot values in intervals of 0.2 in z^niodel- Only simulated spectra for which P^2 > 1% are taken 
into account (4919 over 5000). The dot-dashed lines delimit the region where (Zphot ~ ^model)/(l + ^modcl) < 10%, while the solid line 
corresponds to the case Zpiiot = Zmodel- The photometric redshift accuracy achieved for the whole sample is shown in the left-hand panel 
while that derived for the subsamples of spectra with /ag^'^' < 0.5 and /agn°' ^ shown in the right-hand panels. 



24 /xm (they are shown in Fig. [IJ . We require the simulated 
spectra to be detected in at least five of the observed bands, 
according to the adopted 5a detection limits shown in Ta- 
ble 2] (from Wada et al. 2008). The corresponding Icr limits 
are used to introduce Gaussian fluctuations on the simu- 
lated fluxes and are adopted as the estimated error on the 
resulting fluxes. When a simulated spectrum is undetected 
at a certain waveband than both its flux and the accompa- 
nying error are set equal to half the 5cr detection limit at 
that waveband. 

For all the three sets of simulations, the resulting fluxes 
are fitted using the reference SED templates presented in 
subsection 12.31 A description of the SED models used to 
build the simulations and of the derived photometric accu- 
racy is provided in the following subsections. 

5.1 Simulation I 

The first set of simulated flux measurements (5000 spectra in 
total) have been obtained from the same SED models used 
to construct the reference SED templates, i.e. TAHQ3 model 
for the starburst and ER95 for the AGN. Each SED compo- 
nent is derived by randomly selecting the values of the SED 
parameters within the intervals specifled in subsection 12.31 
The derived photometric redshifts, Zphot, versus the 
simulated redshifts, Zmodci, are shown in the left-hand panel 



of Fig. [7] where the filled circles with error bars represent 
the mean value and the Icr dispersion obtained by fitting a 
Gaussiar0 to the histogram of the Zphot values in intervals 
of 0.2 in ^modci- Only simulated spectra with P^i > 1% have 
been taken into account (they are represented by the small 
dots in the same figure). On average, the accuracy achiev- 
able on (1 -\- z) is close to 10% (a limit indicated by the 
dot-dashed lines in the same figure) below « ~ 2 and above 
z ~ 3, while in the redshift interval 1.8 z ^ 2.8 the dis- 
persion on the photo-z estimates is particularly high. This 
behavior is mainly due to a significant number of points ly- 
ing below the 10%-accuracy region. Interestingly, we found 
that these data correspond to spectra with a dominant AGN 
component, and therefore with almost featureless power-law 
shape. To prove this we have splitted the whole set of sim- 
ulated spectra into two subsamples according to the input 
value of the AGN fraction, i.e. /Xgn ' < 0.5 and /agn ' > 0.5 
and calculated the photometric accuracy for each subsample. 
The results are presented in the right-hand panels of Fig. [T] 
Power-law like spectra have photometric redshifts which are 
systematically and significantly lower than the input red- 
shifts. Once they are removed from the sample the achieved 



The Gaussian fit has been performed using the SKY_STATS IDL 
routine. 
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Figure 8. Photometric redshift accuracy with the AKARI ^¥^7 Deep Survey derived from the set of simulated spectra generated using 
the models of Efstathiou et al. (2000, 2003) and of Efstathiou & Rowan-Robinson (1995) for the star burst /cirrus and AGN infrared 
emissions, respectively. The meaning of the symbols and of the lines is the same as in Fig. [7] For 4859 simulated spectra (out of 5000) 
the fit has provided P^2 > 1%. 



accuracy on (1 + 2:) gets close to 10% over the whole range 
of redshifts probed by the simulation. 

5.2 Simulation II 

The second set of 5000 simulated spectra has been 
constructed using the model of Efstathiou et al. (2000, 
hereafter ERSOO) for the starburst emission while the 
AGN component has been modelled using the prescription 
of ER95, as done before. We have also included a third 
SED component describing the emission from the diffuse 
interstellar dust heated by old stars and/or quiescent star 
formation (infrared cirrus), for which we have followed the 
model of Efstathiou & Rowan-Robinson (2003, hereafter 
ER03). 

According to ERSOO stars form primarily within opti- 
cally thick giant molecular clouds (GMCs) and their model 
provides prescriptions for the evolution of the GMCs ow- 
ing to the ionization-induced expansion of the HII regions. 
The evolution of the stellar population within the GMC 
is accounted for by using the Bruzual & Chariot (1993) 
stellar population synthesis models with a Salpeter Initial 
Mass Function (IMF) and a star mass range of 0.1-125 M©. 
The absorption/emission properties of the dust within the 
GMCs are derived according to the model of Siebenmor- 
gen & Krugel (1992) which assume three different popula- 
tions of dust grains: large grains, small graphite particles 
and PAHs. Starburst galaxies are treated as an ensemble 



of GMCs at different evolutionary stages. The star forma- 
tion rate is modeled with an exponential functional form 
M oc e"'''^^^. Here we set the e-folding time of the star- 
formation rate, tsb, to 20 Myr while the age of the HII 
region phase, AgesB, and the initial optical depth of GMCs, 
To, are treated as free parameters, with AgesB =0 — 72 Myr 
and TO = 50 - 200. 

The infrared cirrus is characterized by lower optical 
depth and lower temperatures of the dust ( ^ 30 K), com- 
pared to what occurs in starburst regions. However, as in 
ERSOO, stars are assumed to have formed within GMCs but 
their evolution is followed well beyond the complete disper- 
sion of the GMCs (i.e. after ~ 72 Myr), when the starlight 
can be absorbed only by the general interstellar medium. 
The input stellar radiation field is computed from the stel- 
lar population synthesis model of Bruzual & Chariot (1993) 
while the grain model of Siebenmorgen & Krugel (1992) is 
used again to derive the absorption/emission properties of 
the dust. The extinction parameter Av regulates the pro- 
portion of the UV to near-IR light that is absorbed by the 
ISM and re-emitted in the far-IR and sub-mm bands. The 
temperature of the interstellar dust is determined by the 
intensity of the radiation field. ER03 characterize this in 
terms of the ratio $ of the bolometric intensity of the radi- 
ation field to that in the stellar radiation field in the solar 
neighborhood. Since the value of $ influences only the far- 
infrared shape of the SED, i.e. the position of the peak of 
the SED (the higher "1>, the hotter the dust and the lower the 
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Figure 9. Photometric redsliift accuracy with the AKARI NEP Deep Survey derived for the set of simulated spectra generated using 
the SWIRE SED templates. The meaning of the symbols and of the lines is the same as in Fig. [7] The simulated spectra for which the 
fit has provided P^2 > 1% are 4060 out of 5000. 



wavelength of the peak), the exact choice of its value is not 
relevant for our purpose and therefore we just set it to 1. Be- 
fore the complete evaporation of the GMC, the non-spherical 
evolution of the molecular cocoon may allow a fraction / of 
the starlight to escape without any dust absorption from the 
GMC. Here we set f = 1 after 3 Myr as in EROS. The age 
of the cirrus component, Agecra,, and the e- folding time of 
the star-formation rate, rem, are treated as free parameters 
and varied consistently with the redshift of the source (i.e. 
always assuring that AgeciR and rem are lower than the 
age of the Universe at the redshift of the source). The visual 
extinction is chosen within the range Ay = — 3. 

Each simulated SED component, i.e. starburst, cirrus 
and AGN, is derived by randomly selecting the values of the 
SED parameters within the intervals mentioned before. The 
fractions of the bolometric luminosity contributed by the 
AGN and by the starburst+cirrus are selected at random 
between and 1. The contribution from the starburst -f cirrus 
is then randomly divided between the two components (i.e. 
starburst and cirrus) ; in this way we guarantee that the final 
sample of simulated spectra include almost 50% of AGN 
dominated SED, as in the previous sample. 

The results on the photometric redshift accuracy are 
shown in Fig. [S] where the meaning of the symbols and of 
the lines is the same as in Fig. [7] The achieved accuracy on 
(1 -|- z) is close to 10% below z ~ 1.8 but decreases signif- 
icantly above z ~ 2, particularly in the interval z ~ 2 — 3. 
Although AGN power-law like spectra are found to be re- 
sponsible for the significant fraction of data points lying be- 
low the 10%-accuracy region (see right-hand panels in the 



same figure), this effect alone does not account for the rel- 
atively low accuracy achieved in the interval 2 2; ^ 3. In 
this redshift range the only valuable redshift indicators are 
the 1.6 fj,m bump and the PAH 3.3 /im emission feature. The 
simulated spectra including a cirrus component are more 
challenging for our reference SED templates to reproduce. 
In fact the latter are meant for fitting relatively young star- 
burst galaxies only, with (eventually) an additional AGN 
component. It is also worth noting that for 2 ^ z ^ 3 the 

1.6 /im bump is shifted to wavelengths 4.8 A ^ 6.4, where 

the AKARI coverage is quite poor (see Fig. [l|. This con- 
tributes to the observed decrease in the photo-z accuracy in 
that redshift interval. 

5.3 Simulation III 

The third simulation is built from the empirical spectra used 
to fit the SED of different types of sources in the Spitzer 
Wide- Area Infrared Extragalactic(SWIRE) survey (see Pol- 
letta et al. 2008 and references therein). The SWIRE tem- 
plate librarjifl contains 25 templates: 3 ellipticals, 7 spirals, 6 
starbursts, 7 AGNs (3 type 1 AGNs, 4 type 2 AGNs), and 2 
composite (starburst -I- AGN). The elliptical, spiral and star- 
burst templates were generated using the the GRASIL code 
(Silva et al. 1998). The ellipticals correspond to three differ- 
ent ages: 2 ,5 and 13 bilion years. The 7 spirals range from 

* http: / / cass.ucsd.edu / SWIRE / mcp / templates / swire_templates. html 
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Figure 10. Upper panel: photometric redshift accuracy, (52/(1 -|- Zmodel); ^ ^ function of the input value of the AGN fraction, /^q^" , 
derived from the three sets of simulations presented here. The filled circles with error bars are the mean value and the Icr dispersion 



provided by the Gaussian fit to the histogram of the &z/{l + ^model) values in intervals of 0.2 in /^q^"^'- 



Only spectra with P 2 > 1% 



have been taken into account. Lower panels: corresponding accuracy achieved on the AGN fraction. Again, the filled circles with error 
bars are the mean value and the la dispersion provided by the Gaussian fit to the histogram of the (/^gn ~ /AGr>f ' ) //agn^ ' values in 
intervals of 0.2 in /ac^ 



early to late types (SO-Sdm). The starburst templates cor- 
respond to the SEDs of NGC6090, NGC6240, M82, Arp220, 
IRAS22491-1808, and IRAS20551-4250. In all of the sphals 
and starburst templates the spectral region between 5 and 
12 /im, where PAH broad emission and silicate absorption 
features are observed, was replaced using observed infrared 
spectra from the PHT-S spectrometer on the ISO and from 
IRS on Spitzer. 

AGN templates corresponding to Seyfert 1.8 and 
Seyfert 2 galaxies were obtained by combining models, 
broadband photometric data (NED) and 750 PHT-S spectra 
of a random sample of 28 Seyfert galaxies. Three other AGN 
templates represent optically selected QSOs with different 
values of infrared/optical flux ratios, derived by combining 
the Sloan Digital Sky Survey (SDSS) quasar composite spec- 
trum and rest-frame IR data of a sample of 35 SDSS/SWIRE 
quasars and then assuming three different IR SEDs. Of the 
remaining AGN templates (type 2 QSOs) one was obtained 
by combining the observed optical/near-IR spectrum of the 
red quasar FIRST J013435. 7-093102 and the rest-frame IR 
data from the quasars in the Palomar-Green sample with 
consistent optical SEDs. The other type 2 QSO template 
corresponds to the model used to fit the SED of a heavily 
obscured type 2 QSO, SWIRE_J104409.95+585224.8 (Pol- 
letta et al. 2006) 

The composite (AGN-I-SB) templates are empirical 
templates created to fit the SEDs of the heavily obscured 
BAL QSO Mrk 231 (Berta 2005) and the Seyfert 2 galaxy 



IRAS 19254-7245 South (Berta et al. 2003). They both con- 
tain a powerful starburst component, responsible for their 
large infrared luminosities, and an AGN component that 
contributes to the mid-IR luminosities. 

Elliptical, spirals, starburst and AGN-I-SB templates 
were randomly selected within the available sample. If the 
SED did not already include an AGN component (i.e. it was 
not one of the two SWIRE composite spectra) than an AGN 
templates was chosen at random from the SWIRE sample 
and added to the SED, with a relative contribution to the 
bolometric luminosity randomly selected from a uniform dis- 
tribution between and 1. 

The results of the photometric redshift accuracy are 
presented in Fig. [5] Again we observe a large scatter in the 
photo-z estimates in the interval 2 ^ z ^ 3 where the main 
infrared features are shifted outside the range of wavelengths 
covered by AKARI and the 1.6 pim feature is not well sam- 
pled because of the lack of coverage around A ~ 5 — 6 /im. 
Despite this, the accuracy achieved on (1 -I- z) is better than 
10% up to 2 = 5. 



5.4 Discussion 

The results of the simulations suggest that for Zmodci ^1-8, 
i.e. when the photo-« estimate is based on the passage of 
the PAH features through the AKARI filters, the achieved 
accuracy on redshift does not depend significantly on the 
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Figure 11. Example of a simulated SED with a dominant AGN 
component (/™Q5\f' > O-^) foi' which the best-fit provides low 
AGN fraction (J^q^ < 0.2). The input SED model is shown 
in the upper panel: the solid line is the total SED and includes 
contributions from a starburst (dashed line) and an AGN (dot- 
dashed), as derived from the models of Takagi et al. (2003) and 
Efsathiou & Rowan-Robinson (1995), respectively. The SED best- 
fit is shown in the lower panel. The meaning of the lines is the 
same as above. 



precise details of the underlying starburst model and the es- 
timated redshift can be considered reliable, irrespective of 
the AGN fraction. At higher redshifts instead, featureless 
AGN-dominated spectra make the recovery of redshift ex- 
tremely challenging from infrared data alone. The effect of 
the AGN infrared emission on the photo-z accuracy is made 
even clearer in the top panels of Fig. [TO] where the quantity 
(zphot - .2modei)/(l + ^modci) IS showu as a fuuctiou of the 
input value of the AGN fraction for the three sets of simu- 
lations presented here. The dispersion on (1 -I- z) increases 
from ^ 10% to ^ 30% as the AGN fraction passes from 

0.1 to 0.5, while for /agn^' ~ 0.5 the photometric estimate 
of the redshift becomes totaly unreliable, at least according 
to the results based on the first two sets of simulations. On 
the other hand, the recovery of the AGN fraction itself be- 
comes very challenging when /agn"' ~ 0.2 as demonstrated 
in the lower panel of the same figure. An example of such a 
case is shown in Fig. 1111 

These results indicates that infrared data alone are a 



valuable tool for redshift estimate but they are not very 
efficient in constraining the AGN fraction, at least when the 
infrared SED is AGN dominated. In such cases the support 
of data at other wavelengths (e.g. optical and near-infrared, 
radio. X-ray) is needed to recover the real nature of the 
infrared emission. 

We note however that AGN power-law like spectra 
can be easily singling out by fitting a line to the observed 
log S — log A diagram [S being the measured flux and A the 
observing wavelength) and analyzing the probability asso- 
ciated to the corresponding , P^\- In Fig. [12] we show 
the photometric redshift accuracy achieved in the three sets 
of simulations after dividing the spectra into those with 
P'^l > 1% (i.e. power-law shape) and those with P'^l < 1%. 
In all cases, the criterion of selection of the SED biased on 
fitting a line to the observed log 5* — log A diagram acts ex- 
actly in the same way as the one relying on the input value 
of the AGN fraction (see right-hand panels of Figs [7] to [9]) . 
Therefore we suggest the following rules as a practical guide- 
line for interpreting photometric redshifts in AKARI NEP 
Deep Survey when using the methods illustrated here: if the 
observed spectrum is poorly described by a power-law than 
the photometric redshift can be considered accurate up to 
2phot ~ 5. When instead the observed spectrum is very close 
to a power-law, it is high probable that the redshift of the 
source has been significantly underestimated and it should 
be rejected. 



6 CONCLUSIONS 

We have tested the reliability of the photometric redshift 
estimates based on infrared photometry alone on a sample 
of 59 galaxies with spectroscopic redshifts 0.5 ^ Zspcc ^1.5 
drawn from the GOODS-N field, combining together in- 
frared data from literature (i.e. Spitzer and ISO) and new 
AKARI 11 and 18 /im photometric data. Our SED fitting 
procedure allows for the contribution to the infrared emis- 
sion from both starburst and AGN, described according to 
the models of Takagi et al. (2003) and Efstathiou & Rowan- 
Robinson (1995), respectively. Three different sets of simula- 
tions derived from both theoretical SED models and empir- 
ical spectra have been used to tested the photo-z accuracy 
achievable with the AKARI NEP Deep Survey. 
The main conclusions are summarized as follows. 

• The SEDs of 42 out of our 59 sources are well fitted (i.e. 
P^2 > 1%) by our starburst+AGN SED reference templates. 
For all the sources in the sample bar 7 the achieved accuracy 
on (1 + z) is close to or better than 10%. Even in the case 
of a "bad" fit (i.e. P^2 < 1%) the mean features in the 
spectrum are recognized and the redshift is still reasonably 
well recovered. Sources at z ^ 2 (three in total) display in 
general a power-law infrared spectrum, whose lack of any 
evident feature make the estimate of redshift quite difficult 
on the basis of infrared data alone. In such cases data at 
longer wavelengths, near or beyond the peak associated to 
the dust emission or at decimetric radio wavelengths, are 
needed to remove the degeneracies in redshifts. 

• Simulations show that the infrared data alone produced 
by the AKARI NEP Deep Survey will provide photomet- 
ric redshifts with a typical accuracy of |2phot — .Zspccl/(1 + 
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Figure 12. Simulated photometric redshift accuracy witli tlie AKARI NEP Deep Survey for power-law like spectra (^^2 > l9o, bottom 
panels) and for spectra which are poorly fitted by a power-law (P^l < t°P panels). The panels refers from left to right to simulations 
based on the Takagi et al. (2003) model for the starbursts, the Efstathiou et al. (2000, 2003) models for the starburst /cirrus and the 
SWIRE sed templates (PoUetta et al. 2008). 



2spec) ~ 10% (Icr) at 2; ^ 2, in agreement with our findings 
for the spectroscopic sample. At higher redshifts the PAH 
features are shifted outside the wavelength range covered by 
AKARI and the 1.6 /nm stellar bump is then exploited as a 
redshift indicator; the accuracy achievable in this case on 
(1 + z) is ~ 10 - 15%, provided that the AGN contribution 
to the infrared emission is subdominant. 

Although our reference SED templates allow for an AGN 
contribution to the infrared emission, the AGN fraction is 
poorly constrained when the SED is AGN dominated and 
therefore does not exhibit any evident features. In these 
cases the support of photometric data at other wavelengths 
(e.g. X-ray, submillimetcr or decimeter radio wavelengths) 
can help to better constrain the relative contributions of 
starburst and AGN to the observed SED. 
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